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q , Abstract 

We predict the energy per baryon in nuclear matter with non-zero frac- 
■ tions of A hyperons. We include Dirac effects on the nucleon as well as the 

Ph ! A and describe how the latter are implemented. We use the nucleon-hyperon 

If} ' meson-exchange potentials from the Jiilich group, the latest as well as an earlier 

version. The dependence of the results on the many-body framework and on 
the nucleon-hyperon interaction model is discussed. 

PACS number(s): 21.65.+f,21.80.+a 
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1 Introduction 

There are important motiva 

presence of hyperons in the interior of neutron stars is reported to soften the equa- 



> 

CNj ' There are important motivations for including strange baryons in nuclear matter. The 

cn ' 

00 

O ■ tion of state (EoS), with the consequence that the predicted neutron star maximum 

. masses become considerably smaller With recent constraints allowing maximum 

masses larger than previously accepted limits [2], accurate microscopic calculations 
which include strangeness (in addition to other effects, such as those originating from 
relativity), become especially important and timely. 

Moreover, as far as terrestrial nuclear physics is concerned, studies of hyperon 
energies in nuclear matter naturally complement our knowledge of hypernuclei (see, 
for instance, Refs. [3j [51 [6j 0] . For example, the EoS of hypermatter is useful in the 
development of generalized mass formulas depending on density and hyperon fraction 
[7J Hj. From the experimental side, additional data are very much needed, especially 
on AA hypernuclei, which would provide information on the A — A interaction. Con- 
cerning single hypernuclei, analyses of data on A binding energies constrain the depth 
of the single-A potential to be 27-30 MeV [8]. The status of S hypernuclei and the 
E-nucleus potential is more controversial (see Ref. [9] and references therein). 

Microscopic calculations of nuclear matter properties including hyperons have 
been reported earlier within the non-relativistic Brueckner-Hartree-Fock framework 
(BHF) (see, for instance, Refs. [TJ, [10] ) , using the Nijmegen [11] nucleon-hyperon (NY) 
meson-exchange potential. Extensive microscopic work on hyperonic nuclear matter, 
again within the non-relativistic BHF framework, has also been published by the 
Barcelona group (see, for instance, Refs. [141 [To] [To ] fTTj). 



1 



It is one purpose of the present work to bring in the new aspect of Dirac ef- 
fects on the A hyperon as well as the nucleon. By "Dirac effects" we mean that the 
single-baryon wavefunction is calculated self-consistently with the appropriate effec- 
tive interaction. The importance of these effects on the nucleonic equation of state 
cannot be overestimated [TB] . They can be seen both as relativistic effects and three- 
body forces originating from nucleon-antinucleon excitations, and provide an essential 
saturation mechanism missing from conventional approaches. 

Our previous calculation [19] of the binding energy of a A impurity in nuclear 
matter showed that Dirac effects on the A hyperon yield a moderate reduction of the 
binding energy. In that calculation, we used the most recent nucleon-hyperon (NY) 
potential reported in Ref. [20J (thereafter referred to as NY05), and observed that 
G-matrix predictions obtained with NY05 are significantly differerent from calcula- 
tions using the previous (energy independent) version of the Jiilich NY potential [13J. 
Therefore, in this work we will use both potentials, for comparison. The Bonn B 
potential [21] is used throughout for the nucleon-nucleon (NN) part. 

Previous calculations of the EoS with hyperons have typically been conducted 
within a non-relativistic framework together with r-space local NY potentials. Our 
approach uses non-local relativistic momentum-space (NN and NY) potentials and a 
relativistic many-body method, and is therefore fundamentally different. 

The baseline work for our EoS's was developed in Ref. [22]. The framework 
described there for isospin-asymmetric nuclear matter has been further expanded 
and adjusted to include other species of baryons in different concentrations. Our 
latest EoSs for the NN sector can be found in Ref. [23] . 

In the next Section, we describe some technical aspects of the calculation. We 
then proceed to present and discuss our results (Section 3). Conclusions and future 
plans are summarized in Section 4. 



2 Description of the calculations 
2.1 General formalism 

For a total density of baryons (nucleons and A's), p, and some specified A fraction, 
1a, the densities of each species are known and are related to their respective Fermi 
momenta by 

/^ = ^(2)(^) 3 , Pa = ^(^) 3 - (1) 
For a nucleon and a A with momenta k^ and k A , the total and relative momenta 



are 



p t l? I MNkA ~ M ^ N I rZ fo\ 
P=k N + k A , k = — — — — = ak A -pk Nl (2) 
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with a and (3 equal to M ^ M and M ff+ M , respectively. Thus, 

k A = f3P+ k, k N = aP- k. (3) 
Clearly, for the case of two nucleons we have 

— * — * 

P = k N + k' N , k= k -^L. (4) 

The single-nucleon and single-A potentials are obtained as 

U N (k N ) = U NA (k N ) + U NN (k N ), (5) 

and 

U A (k A ) = U AN (k A ), (6) 

i.e., the AA interaction is neglected. In the equations above, the various terms, Ub 1 b 2 , 
are the contributions to the potential of baryon B\ from its interation with all baryons 
of type B 2 . They are given by 

U NA (k N )= E if x 1 1 wo J + 4-\\ G T A' S ' J (Hk N , k A ),P(k N , k A ))d"k A , (7) 

^ra 7 2tjV + 1 )[2s N + 1 JO 



-^ S J (2t N + l)(2s N + 1) 

(2r+l)(2J+l) 
, LSJ (2t N + l)(2s N + l) Jo 



u NN (k N )= E if x J x^ /*" QS£ ,g,J (fc(^. (8) 

rrTi t 2i/v + 1 )(2sn + 1 JO 



and 



^(4)= e if x y + a r G T A ^ j (k(k N , k A ),p(k N , k A ))d 3 k N , 



(2*a + 1)(2sa + 1) 



where the channel isospin T can be or 1 for the NN case and is equal to 1/2 for the 
NA case, and Si,U (i = N,A) are the spin and isospin of the nucleon or A. 
Notice that 

P- « ^ (10) 

U AN PN 

an approximation often used in mean-field approaches. 

The average potential energies of nucleons and A's are determined from 

< Un >= —7^4 r u N& N )dk? N , (11) 
Pn (27r) d 2 Jo 

and 

111 /"^ A 

< Ua >= 7T7^ 2 9 / ' MkA)dk\, (12) 
Pa (^7t; 2 jo 
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where the factors of 4 and 2 in Eqs. (11) and Eq. (12), respectively, account for 
protons and neutrons in both spin states or A's in both spin states. 
Finally the average potential energy per baryon is obtained as 

<u>= Pn<Un>+P,<Ua>^ (13) 

Ptot 

from which, together with a similar expression for the kinetic energy, one obtains the 
average energy per baryon. 

The NA G-matrix is obtained from the Bethe-Goldstone equation 

< NA\G NA (E )\NA > = <NA\V\NA> (14) 

+ J2 < NA\V\NY >— 0— <NY\G NA (E )\NA>, 
y=A,s Hjq-h, 

where E Q and E are the starting energy and the energy of the intermediate NY state, 
respectively, and V is an energy-independent NY potential. 

For two particles with masses and M A and Fermi momenta kp and kp, Pauli 
blocking requires 

n(hP\-\ l \PP + k\ > kp and \aP - k\ > kp 1 n * 

WjPj "\0 otherwise. (15j 

The above condition implies the restriction 

<MN P y + k 2_( k Ny ( A^p )2 + k 2 _ (fc A)2 



2Pk^ >cosu> ^ p , Ubj 



where 9 is the angle between the total (P) and the relative (k) momenta of the two 
particles, and M = M A + M^. Angle-averaging is then applied in the usual way. 

In the present calculation we consider a non-vanishing density of A's but do not 
allow for the presence of real E's in the medium (although both A and S are included 
in the coupled-channel calculation of the NY G-matrix, see Eq. (14)). Essentially we 
are considering a scenario where a small fraction of nucleons is replaced with A's, 
as could be accomplished by an experiment aimed at producing multi-A hypernuclei. 
Multistrange systems, such as those produced in heavy-ion collisions, may of course 
contain other hyperons on the outset. For small A densities, though, as those we 
consider here, the cascade (5) and the S hyperon are expected to decay quickly 
through the strong processes + H —>■ A + A and N + S — ► + A. Under these 
conditions, a mixture of nucleons and A's can be considered "metastable" , in the 
sense of being equilibrated over a time scale which is long relative to strong processes 
[7]. (In fact, the strong reactions mentioned above would have to be Pauli blocked in 
order to produce a metastable multistrange system [24].) 

We neglect the YY' interaction, as very little is known about it. Furthermore, 
non-local momentum-space YY potentials, appropriate for our DBHF framework, are 
not available at this time. For these reasons, we keep the A concentration relatively 
low. 
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2.2 Dirac effects on the NA potential 

The relation between the non-relativistic Brueckner approach and the relativistic 
framework (known as Dirac-Brueckner-Hartree-Fock, DBHF) has been discussed for 
a long time. Already in Ref. it was shown how relativistic effects tie in with virtual 
excitations of pair terms. Lately, these concepts have been revisited [26J with similar 
conclusions. In short, the Dirac effect on the EoS of nucleonic matter is an essential 
saturating, and strongly density dependent, mechanism, which effectively accounts for 
the class of three-body forces originating from virtual nucleon-antinucleon excitations. 
When hyperon degrees of freedom are included, for reasons of consistency, those 
should then be subjected to the same correction. 

We have incorporated DBHF effects in the A matter calculation, which amounts 
to involving the A single-particle Dirac wave function in the self-consistent calculation 
through the A effective mass, M£. Similarly to what is done for nucleons in Ref. [2Tj . 
we fit the single-particle energy for A's using the ansatz 



e A (p) = v /(Ml) 2 +p 2 + ^, (17) 

where M A = Ma + Ug, and Ug and Uy are the scalar and vector potentials of the 
A baryon, which we assume to be momentum independent. Note the single-particle 
energy is the sum of the single-particle potential, for which we use the two-parameter 
ansatz, 

M* 

U A (p) = Ug + U$, (18) 

'{M* A ) 2 +p 2 



and the expectation value of the free energy operator in the Dirac equation [21] . 
Because of the small hyperon concentrations considered here, the parameters of the 
single nucleon potential are taken from a separate calculation of symmetric nuclear 
matter J22J 122] performed at densities corresponding to the nucleon densities, 
in Eq. (1). This amounts to saying that, for low hyperon concentrations, the main 
correction to the single nucleon potential comes from the presence of the first term 
in Eq.(5). The parameters of the single A potential are fitted, self-consistently with 
the NA G-matrix, at momenta close to its Fermi surface, specifically k\ = 0.8kp and 

fc 2 = fc£. 

A problem with the Jiilich NY potential in conjunction with DBHF calculations 
is the use of the pseudoscalar coupling for the interactions of pseudoscalar mesons 
(pions and kaons) with nucleons and hyperons. For the reasons mentioned above 
(that is, the close relationship between Dirac effects and "Z-diagram" contributions), 
this relativistic correction is known to become unreasonably large when applied to 
a vertex involving pseudoscalar coupling. On the other hand, the gradient (pseu- 
dovector) coupling (also supported by chiral symmetry arguments) largely suppresses 
antiparticle contributions. To resolve this problem, one can make use of the on-shell 
equivalence between the pseudoscalar and the pseudovector coupling, which amounts 
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to relating the coupling constants as follows: 

9ps fpv ) (19) 



where g ps denotes the pseudoscalar coupling constant and f pv the pseudovector one; 
m ps , Mi, and Mj are the masses of the pseudoscalar meson and the two baryons 
involved in the vertex, respectively. This procedure can be made plausible by writing 
down the appropriate one-boson-exchange amplitudes and observing that, redefining 
the coupling constants as above, we have (see Ref. [21] for the two-nucleon case) 

V pv = V ps + (20) 

where the ellipsis stands for off-shell contributions. Thus, the pseudoscalar coupling 
can be interpreted as pseudovector coupling where the off-shell terms are ignored. 
This is what we apply in our DBHF calculations. 

In the coupled channel calculation, evaluation of Gna involves the transition po- 
tentials Vjva^ata, VjvA<->-iV£) an d Vnz^nz (ah with total channel isospin equal to 1/2), 
plus the corresponding exchange diagrams. Because in the present scenario the E 
hyperon is not given an effective mass, Dirac effects are applied only in Vata->jva- A 
diagram where not all of the baryon lines are Dirac-modified may yield a Dirac effect 
that is artificially skewed. 

Finally, a comment is in place concerning meson propagators. In standard DBHF 
calculations [21J, the so-called Thompson equation (a relativistic three-dimensional 
reduction of the Bethe-Salpeter equation) is used for two-baryon scattering. In the 
Thompson formalism, static propagators are employed for meson exchange, i.e. 

(21) 



where m a denotes the mass of the exchanged meson and q, <f are the baryon momenta 
in their center-of-mass frame before and after scattering. The Julich NY potentials 
are based upon time-ordered perturbation theory [T2] and use a meson propagator 
given by 

—( f^f \ ( 22 ) 

LO a (z - Ei — Ej - u a ) 



with u a = \jvp? a + (<f — q) 2 \ Ei = y M? + <f 2 and Ej = J Mj + q 2 are baryon energies 
and z is the starting energy of the two-baryon system. In order to eliminate the energy 
dependence, Reuber et al. [13] replaced the original z with 

z=^(M 1 + M 2 + M 3 + M 4 ), (23) 

where the Mj's denote the baryon masses of the four legs in the one-meson exchange 
diagram. In any case, the Julich meson propagator involves the baryon masses. Re- 
placement of these free-space masses with in-medium values would create medium 
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effects on meson propagation which we do not wish to include in our nuclear matter 
calculations. The reason for keeping free-space masses in the meson propagator is 
twofold. First, standard DBHF calculations do not include medium effects on meson 
propagation as they typically use Eq. (19), which does not depend on baryon masses. 
Second, medium effects on meson propagation constitute a separate class of effects 
that we are not concerned with in the present context and are typically not perceived 
as part of the DBHF approach. 

Having taken the steps described above, we proceed to the DBHF calculations. 
The effect of such corrections on the binding energy of a A impurity in otherwise pure 
nuclear matter was reported in Ref . [IH] . 

3 The equation of state of nucleon and A matter 

As stated in the Introduction, we will be showing results for both the NY05 potential 
[2"U] and the previous version of the Julich NY potential, NY94, specifically the model 
referred to as A in Ref. |13j . 

In Fig. 1 we show the energy per particle as a function of density for different A 
concentrations as obtained from DBHF calculations along with the NY94 potential. 
As more nucleons are replaced with A's, generally less binding energy per particle is 
generated. This is due to the weaker nature of the NA interaction relative to the 
NN one. Furthermore, the A Fermi momentum grows rather quickly with p\, since 
only two A's can occupy each state, rather than four, which implies a fast rise of the 
hyperon kinetic energy. (Although, for small hyperon densities, there is at first some 
reduction of the kinetic energy, due to the fact that A's have larger mass and can 
occupy lower energy levels.) We notice that the saturation density remains essentially 
unchanged with increasing hyperon concentrations. As density grows, though, larger 
A concentrations start to yield increased attraction. One must keep in mind that, 
especially at the higher densities, the NN interaction become less and less attractive 
due to medium effects, in particular repulsive Dirac effects. Thus, removing nucleons 
from the system can actually amount to increased binding. 

Moving now to Fig. 2, where the NY05 potential is adopted, we see a very different 
scenario. We recall that the NY05 model is considerably more attractive, yielding 
about 50 MeV for the A binding energy [271 EH] whereas a value close to 30 MeV 
was found with NY94 [13J. Naturally, we expect these differences to reflect onto 
the respective EoS predictions. Here, increasing the hyperon population yields more 
binding, a trend opposite to the one seen in the previous figure. Again, the final 
balance is the result of a combination of effects. The fact that the NN component is 
repulsive at the higher densities, together with the more attractive nature of the NY05 
potential, determines here a net increase in attraction with decreasing nucleon density. 
However, the additional binding becomes smaller and smaller with increasing hyperon 
concentration, indicating that, at sufficiently large A densities, the net balance may 
turn repulsive. Notice also that in Fig. 2 the minimum moves towards higher densities, 
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signifying that baryon pairs favor a smaller interparticle distance. 

To summarize, one must keep in mind that the energy/particle is the result of 
a delicate balance of both the kinetic energies and the potential energies of the two 
baryon species, being weighed by the respective densities. Thus, although the NN 
interaction is generally more attractive than the NA one, the net effect of replacing 
nucleons with A's will depend sensitively on the nature of the NY interaction that's 
being put into the system, as well as the "stiffness" of the original, nucleonic, EoS. 

This is further confirmed in Fig. 3, where we show predictions from conventional 
BHF calculations (i.e., no "Dirac" effects). As in the previous figures, the red curve 
is the nucleonic EoS. The solid(dashed) curves are the predictions with NY94(NY05). 
Clearly, the effects are opposite depending on the NY interaction. Qualitatively, the 
trend seen in each group of curves (solid or dash) is approximately consistent with 
the one observed previously in the corresponding (i. e. same NY potential) DBHF 
predictions. However, the "cross over" of the curves seen in Fig. 1 at about twice 
saturation density is due to a large extent to the more repulsive nature of the nucleonic 
EoS in the DBHF calculation (see comments above). 

In conclusion, the effect of hyperons on the EoS is strongly dependent upon the 
baseline (nucleonic) EoS as well as the NY potential model. With regard to the 
first issue, we stress the importance of starting from a realistic EoS, such as the one 
predicted in DBHF calculations or obtained with the inclusion of three-body forces. 
The second observation confirms the conclusions of Ref. [25] concerning the large 
uncertainties originating from the bare NY potentials. Unfortunately, the existing 
data do not set sufficient constrains on the potentials, as demonstrated by the fact 
that different potentials can fit the available scattering data equally accurately but 
produce very different scattering lengths [28J. 

4 Conclusions 

We have reported on Dirac-Brueckner-Hartree-Fock predictions of the energy per 
particle in symmetric nuclear matter as a function of the total baryon density and 
(moderate) A concentrations. Dirac effects are included on both the NN and the 
NA potentials. Our DBHF scheme, which requires the use of relativistic momentum- 
space nucleon-baryon potentials, represents a fundamentally different paradigm as 
compared to existing calculations. 

Ultimately, the actual fraction of hyperons present in star matter must be deter- 
mined by the equations of (3 stability and charge neutrality for highly asymmetric 
matter containing neutrons, protons, hyperons, and leptons. Unlike the straightfor- 
ward calculation of lepton fractions in /3-stable isospin-asymmetric nucleon matter, 
the case of nucleons and hyperons require a lot more effort due to their strongly inter- 
acting nature. One needs to know the chemical potentials of each species (neutrons, 
protons, various hyperons) ideally at any densities and baryon concentrations. Those 
may then have to be fit with analytic functions (not a trivial task, given the num- 
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Figure 1: (Color online) Energy /particle as a function of density in symmetric nuclear 
matter for various A concentrations Y\. Predictions obtained from DBHF calculations 
with the NY94 potential. 
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Figure 2: (Color online) Energy/particle as a function of density in symmetric nuclear 
matter for various A concentrations Predictions obtained from DBHF calculations 
with the NY05 potential. 
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Figure 3: (Color online) Energy /particle as a function of density in symmetric nuclear 
matter for various A concentrations. Predictions obtained from BHF calculations. 
The solid(dashed) lines are obtained with the NY94(NY05) potentials. 



ber of dependences), so that they can be inserted into the appropriate equations to 
determine the equilibrium fractions of each baryon type. This is deferred to another 
paper. 

We observed that the predicted energy/particle is extremely sensitive to the chosen 
NY interaction. The large uncertainties due to the model dependence discussed in this 
paper are likely to impact any conclusions on the properties of strange neutron stars, 
which therefore must be interpreted with caution. These include considerations of 
deconfmement and possible transition from hadronic to quark matter, which depend 
sensitively on the equation of state in the hadronic phase. 
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